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A model is introduced describing the interplay between superconductivity and spin-ordering. 
It is characterized by on-site repulsive electron-electron interactions, causing antiferromagnetism, 
and nearest-neighbor attractive interactions, giving rise to d-wave superconductivity. Due to a 
special choice for the lattice, this model has a strong-coupling limit where the superconductivity 
can be described by a bosonic theory, similar to the strongly coupled negative U Hubbard model. 
This limit is analyzed in the present paper. A rich mean-field phase diagram is found and the 
leading quantum corrections to the mean-field results are calculated. The first-order line between 
the antiferromagnetic- and the superconducting phase is found to terminate at a tricritical point, 
where two second-order lines originate. At these lines, the system undergoes a transition to- and 
from a phase exhibiting both antiferromagnetic order and superconductivity. At finite temperatures 
above the spin-disordering line, quantum-critical behavior is found. For specific values of the model 
parameters, it is possible to obtain SO (5) symmetry involving the spin- and the phase-sector at 
the tricritical point. Although this symmetry is explicitly broken by the projection to the lower 
Hubbard band, it survives on the mean-field level, and modes related to a spontaneously broken 
SO (5) symmetry are present on the level of the random phase approximation in the superconducting 
phase. 
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I. INTRODUCTION 

Both for empirical- and historical reasons, research on 
superconductivity tends to be preoccupied with the weak 
coupling limit. From a more general perspective, BCS 
theory as well as Gorkov-Migdal-Eliashbcrg theory corre- 
spond with a special case which in a sense is pathological. 
The emphasis is completely on the amplitude of the or- 
der parameter while fundamentally superconductivity is 
about breaking of gauge symmetry, associated with the 
phase sector. The work of Schmitt-Rink and Nozieres 
revealed that the BCS theory for a s-wave supercon- 
ductor can be smoothly continued to the strong coupling 
limit. It is generally recognized that it is far easier to 
understand the vacuum structure of such a superconduc- 
tor in strong coupling. Amplitude fluctuations can be 
regarded as highly massive excitations and all what re- 
mains is the phase sector described in terms of hard-core 
bosons, or alternatively in terms of pseudo-spin models. 

In the context of high T c superconductivity one en- 
counters a far more complex physics. Abundant evidence 
is available for a d-wave superconducting order parame- 
ter. This is usually discussed in terms of weak-coupling 
theory with its d-wave nodal fermions while the more so- 
phisticated approaches start from this limit, attempting 
to penetrate the intermediate coupling regime using self- 
consistent perturbation theory Q . The obvious problem 
is that the coherence length is rather short Q. At the 
same time, an interesting case has been presented claim- 
ing that much of the thermodynamics can be understood 
from phase-dynamics alone B, completely disregarding 



amplitude fluctuations. It would therefore be useful to 
study strong coupling theories for d-wave superconduc- 
tors. 

An even better reason to pursue a strong coupling per- 
spective is the growing evidence for the presence of well 
developed antiferromagnetism coexisting with the super- 
conductivity. Traditionally, this was approached within, 
again, an implicitly weakly coupled perspective. The 
magnetic fluctuations as seen in NMR and neutron scat- 
tering were believed to be due to the proximity to an 
amplitude driven spin density wave transition Re- 
cently, this perspective has been drastically changed due 
to the observation of strong static antiferromagnetic or- 
der associated with the stripe phases in the La,2CuO/± 
system || . In the Nd doped samples where this order is 
strongest the magnitude of the Neel order parameter can 
be as large as 0.3 [is 0, while 0.1 fis has been claimed 
in 'pristine' Lai.ssSro.^CuO^ ||. It appears that this 
antiferromagnetic order is in competition or even coexist- 
ing with the superconducting order [^|j9| . Given that the 
stripe antiferromagnet should be strongly renormalizcd 
downward due to transversal quantum spin fluctuations 
fiof the stripe antiferromagnet has to be strongly cou- 
pled. Given the strong similarities between the static 
order and the incommensurate spin fluctuations which 
seem to be generic for all cuprate superconductors in the 
underdoped regime, a strong coupling perspective on the 
antiferromagnetism should be closer to the truth even if 
static order is not present, at least as long as the doping 
is not too large. 

Recently several theoretical attempts have been under- 



1 



taken to shed light on this problem of strongly coupled 
superconductivity and antiferromagnetism. The simplest 
theory of this kind is Zhang's £0(5) theory, where super- 
conductivity and antiferromagnetism are 'unified' within 



a single larger symmetry 1 1 1 . Given that no such symme- 
try is manifestly present at the ultraviolet of the prob- 
lem, this might well be misleading and one would like 
to have a more general framework in which this (near) 
50(5) symmetry appears as a special case. The manifest 
symmetry of the problem is U(l) X SU (2) (superconduct- 
ing phase- and spin, respectively). The structure of the 
long wavelength effective theory based on this symmetry 
principle has been analyzed recently by one of the au- 
thors [|l2| , including the charge order associated with the 
stripe phase. These approaches are only truly meaningful 
at long wavelength and a more complete understanding 
is in high demand. In fact, the only reasonably complete 
theory is the one by Vojta and Sachdev , based on the 
large iV/small £ saddle point of the Spl(2N) t— J model. 
However, in this large N limit the antiferromagnetism is 
in the strongly quantum disordered regime, and is there- 
fore at best dual to the renormalized classical Neel order 
of the stripe phases. 

Here we will present an exceedingly simple toy model 
which seems nevertheless to catch much of the physics 
discussed in the above. It is similar in spirit to the 
lattice-boson description of superconductivity and mag- 
netism discussed in Ref's Q and ]l5| ]. The pursuit is 
to construct a model which at the same time describes 
localized magnetism and local pairing superconductivity. 
The magnetism is undoubtedly related to strong, Hub- 
bard U type on-site repulsions. This prohibits for obvi- 
ous reasons on-site paring. The next microscopic length 
scale available on the lattice is the lattice constant itself: 
the pairs causing the superconductivity live on the links 
of the lattice |lj|. If such a link-pair is occupied, the 
sites connected by this link are both occupied by a single 
electron. In the presence of on-site repulsions these elec- 
trons will tend to turn into a spin system. The number 
fluctuations implied by the superconducting phase order 
correspond with such an occupied link-pair becoming un- 
occupied, causing at the same time a dilution of the spin 
system. 

On the square lattice a subtlety keeps a theory with 
these link pairs as building blocks from being simple. Dif- 
ferent from the large N limit with its spin-Peierls order 
Jk| , the link pairs cause both conceptual problems in de- 
scribing the state at half-filling as well as serious technical 
problems. As will be discussed in Section II, a consistent 
formulation requires local constraints to be added to the 
theory in order to exclude tilings of the lattice charac- 
terized by multiple occupancies on the sites. This is not 
necessarily fatal: the theory is bosonic and it might well 
be that Jastrov projections cure the problem. A central 
result of this paper is our discovery of a different lattice 
where these likely non-essential 'correlation' problems are 



absent: the 1/5 depleted lattice shown in Fig. 1. The 
link-pairs live on the long bonds, while the short bonds 
only carry spin-spin interactions. As will be further dis- 
cussed, this model is characterized by an unproblematic 
classical (in fact, large d) limit. This allows us to derive 
in a controlled way a complete semi-classical description. 

As discussed in Section III, we find a surprisingly rich 
phase diagram on the classical level containing all phases, 
which have been up to now suggested in this context, 
including the large N spin-quantum paramagnets. Per- 
turbing around this classical limit, we address the struc- 
ture of the semiclassical theory including the universality 
classes at the various phase transitions (Section IV). By 
fine tuning parameters, we find lines in the phase dia- 
gram where the £0(5) symmetry is approached. How- 
ever, even at the most symmetric point £0(5) is not 
reached: as we will show, the theory becomes £0(5) sym- 
metric on the classical level but the quantum corrections 
destroy this symmetry again. As was already pointed out 
in the context of the £0(5) symmetric ladders, fine tun- 
ing of the on-site repulsions is required to stabilize the 
full symmetry (Sections V and VI). 



II. THE MODEL 

A. Correlated superconductivity 

For the strong-coupling description we are aiming at, 
the microscopic building blocks are electron link-pairs, 
created by the operators 
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where S is a lattice unit-vector, while i labels the sites. 
Such a link-pair is the typical microscopic object in a 
strong-coupling theory of d-wave superconductivity and 
the smallest electron-pair that can support spin degrees 
of freedom. Two serious problems arise when trying to 
construct a model from these operators, one technical 
and one conceptual. The technical problem is related 
to the spatial structure of the link-pairs, which intro- 
duces correlations between pairs centered on different 
bonds. These correlations show up in the commutation 
relations of the link operators. Operators along differ- 
ent bonds do not commute if their links share a common 
site. As a result, the dimension of the link-operator al- 
gebra grows with the system size. This makes a simple 
pseudo-spin description of the charge-sector impossible 
and not much seems to have been gained by going to the 
strong-coupling limit. 

This problem can be avoided by assuming that one 
can somehow keep track of which electrons belong to a 
particular pair (this can be ambiguous, for instance in 
the case of four electrons sitting in a square). If this 
is possible, the link-pairs can be described by hard-core 
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boson-operators, satisfying bfV^&fxf 2 ^ = for S 7^ 5'. 
Link-bosons on different bonds always commute, remov- 
ing the problem of the infinite-dimensional link-algebra. 
The correlation effects then show up in a different way, 
however. The hard-core link-bosons are spinful gener- 
alizations of the quantum dimers (IT]] pl|. It is well 
known that even the classical theory of the dimers is a 
complex combinatorics problem, which was solved for the 
case of half-filling but not for general densities. This 
problem seems unavoidable when one tries to construct a 
strong-coupling theory for electron-pairs with one or the 
other real space internal structure on the square lattice. 

The conceptual problem is related to the fact that our 
link-pairs carry spin. It concerns the state at half-filling. 
On the square lattice, there are many ways in which the 
link-pairs can be distributed over the lattice to obtain 
complete covering. Since the half-filling state is a pure 
spin-system, this charge degree of freedom is superfluous. 
The link-pair model at half-filling therefore suffers from 
a large degeneracy. 

In the large- N t-J model studied by Vojta and Sachdev 
p3| , link-pairing arises as a result of nearest-neighbor 
spin-singlet formation, and the pairs are in this case spin- 
zero dimers. As a result, different link-pair configurations 
at half-filling correspond to different distributions of the 
singlet spin-bonds over the square lattice. These configu- 
rations are therefore physically distinct. The spin-Peierls 
order which is present at half-filling singles out a partic- 
ular link-pair configuration, breaking the degeneracy. 

For a large S type antiferromagnet, however, the spin- 
sector cannot be used to break the degeneracy associated 
with half-filling. Let us therefore consider a model where 
link-pairing arises as a result of charge-charge interac- 
tions. In this case, link-pairs can have both a singlet 
and a triplet spin-component, allowing for the construc- 
tion of a half- filling antiferromagnet. Consider a nearest- 
neighbor attractive interaction V, an on-site repulsive in- 
teraction U and a longer-range repulsive interaction U' , 
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hopping processes , (2) 



where S runs over all lattice unit-vectors. The attrac- 
tive interaction V promotes link-pairing, while the longer 
range repulsive interaction U' is needed to counteract 
phase-separation in the strong-coupling limit. 

At small electron densities, the strong-coupling limit 
of the above model describes a dilute gas of electron 
link-pairs. Near half-filling, it describes a dilute gas of 
hole link-pairs, moving through a spin background. Tak- 
ing hole-pairs and spins, instead of electron-pairs, as the 
elementary building blocks in the strong-coupling limit 
near half-filling, the large degeneracy in the description 



is avoided. Such a perspective is not entirely satisfactory, 
however, since the spin-sector is in this case represented 
in a first-quantized form. 

The technical problems, related to the spatial correla- 
tions between the link-pairs, of course remain also for this 
model. These correlations become important at finite 
densities away from zero- or half-filling, severely compli- 
cating the strong-coupling analysis of this model. More- 
over, the short-range attractive and long-range repulsive 
interactions will give rise to charge ordering phenom- 
ena at intermediate densities, further complicating the 
physics. 




FIG. 1. The ^-depleted lattice. Dotted lines connect near- 
est-neighbor horizontal and vertical bonds. 



B. Depleted lattice 

The complex spatial correlations between link pairs 
and the tendency towards charge-ordering at intermedi- 
ate densities as discussed in the previous subsection can 
be avoided by formulating the model not on the square 
lattice, but on the 1/5-depleted lattice, shown in Figj^. 
We arrive at this lattice by expanding the sites of a 
square lattice to form tilted squares. Along the bonds 
of the original square lattice, attractive charge-charge 
interactions are assumed, while on-site repulsive inter- 
actions are introduced to promote antiferromagnetism. 
The electron-Hamiltonian of such a model reads 



H 



n- 



t t i.S i 

Vn{ n 2 



T n i| 



n 2 | n 2 ^ 



-hopping processes , 



(3) 



where the index i labels the square plaquettes, while (i, 5) 
denotes the four bonds extending from these plaquettes. 
The two sites connected by each long bond are numbered 
1 and 2 from left to right and from bottom to top. The 
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hopping processes can include hopping along the long 
and the short bonds, as well as longer-range hopping 
across the square or the octagonal plaquettes. In the 
large V, large U limit, the above model reduces to one 
describing the physics of spinful link-pairs, which reside 
on the long bonds of the 1/5-depleted lattice. Note that 
the spatial correlations between these pairs are the same 
as between point-particles on a square lattice. Since the 
link-pairs on different long bonds do not share a common 
site, the algebra of the link-pairs on different bonds de- 
couples and a pseudo-spin type description of the charge 
sector becomes possible. Admittedly, this amounts to a 
rather radical simplification as compared to the square 
lattice link-pair problem. However, the long wavelength 
physics we will derive for the depleted lattice might be 
of a greater generality because of the universality princi- 
ple. In fact, we suspect that the complexities discussed in 
the previous subsection will add only tendencies towards 
charge ordering which can be to some extent discussed 
separately. 



C. Pair- hopping and spin-spin interactions 
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where the same notation has been used as in Eq. 
A projection operator Pa = (1 — n^n'^fl — n|*nfjj 
has been included in the definition of the link-operators 
L^g" 2 *, Eq. (|l|). This enforces the constraint of no dou- 
ble occupancy, which is a result of the large U limit in 
Eq.®. 




Since the Hamiltonian Eq. (|3j) should be viewed as a 
toy-model, there is no reason to explicitly derive the 
strong-coupling description by starting from this Hamil- 
tonian and integrating out the states with unpaired 
electrons. Instead, we simply formulate another toy- 
model, which describes generic features of the dynamics 
of bound link-pairs on the 1/5-depleted lattice. We in- 
clude the minimal number of processes needed to capture 
the physics of such a system, making sure that the inter- 
actions are consistent with the symmetries of the lattice. 

An antiferromagnetic spin-spin interaction J is as- 
sumed along the long bonds and a ferromagnetic inter- 
action Jp along the short bonds. This choice allows for 
an extension of the model to higher dimensions without 
introducing frustration into the spin system, making it 
possible to reach the d — > 00 limit and check the mean- 
field results there. For Jp ^S> J, the half-filled system 
becomes equivalent to an S = 2 antiferromagnct on a 
square lattice (or S — d on a d-dimensional hypercubic 
lattice). This property will be used to obtain an estimate 
of the quantum-corrections to the saddle-point results ob- 
tained in the next section. 

A sublattice and an inter-sublattice hopping process 
are introduced, with amplitudes t\ and t 2 . Both pro- 
cesses move a pair from a horizontal (vertical) bond to 
a nearest-neighbor vertical (horizontal) bond. The t\ 
process respects the spin-ordering, keeping the electrons 
which form the pair on their original sublattice. The t 2 
process moves the electrons from one sublattice to an- 
other, thereby frustrating Neel order. 

Including a chemical potential /i, we arrive at the 
Hamiltonian 




FIG. 2. Hopping processes and spin-spin interactions in- 
cluded in the model. 

The Hilbert space on one long bond is spanned by five 
states: unoccupied (V), spin-singlet (A) and spin-triplet 
(1,0,-1). The operators acting on this space are 5x5 
matrices. Introducing the notation 



the pair creation operators can be written as: 
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These operators are the equivalent of the pseudo-spins 
which appear in the strong-coupling negative U Hubbard 
model pQ ]. The operators G Q y, Gy a and \(n a — ny) 
form an S — \ spin-algebra (a = 1,0,-1, A). Pseudo- 
spins with a different spin-index a do not commute. In 
section the constraint of no double occupancy is aban- 
doned to allow for the construction of an SO(5) symmet- 
ric version of this model. The operators (pi) then become 
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S = 1 pseudo-spins and operators with a different index 
a do commute in this case. 

It is convenient to introduce the total spin and the Neel 
moment of a link-pair 

5*1,(5 = su,s + s*2i,a ; Si ; s — su t g — s*2i,<5 , (7) 
which are given by 

O — On — Lr-1-1 ; 

5+ = \/2(Gio + Go-i) , 

D — — UAO — U-OA ; 

5+ = >/2(Gia-Ga-i), (8) 

satisfying 50(4) commutation relations. After absorb- 
ing a factor (— iy*+ l v sign(S) into the singlet-state l-A^j), 

which induces a staggering of 5 and Gav, the Hamilto- 
nian takes the form 



« = E E 



t 5 X = ±6^ 

5 2 = ±<Sy 



ce=l,0,-l 

+(ti-t 2 )(G^G^ + h.c.) 



EE 
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1 * 

1 — n. 



(9) 



where % = (— l)**+ l y is the AF staggering factor. Note 

that it cannot be absorbed into Si t s, since (i,5 x ) and 
(i + 1, — S x ) label the same bond. 

III. MEAN-FIELD ANALYSIS 

A variational Hartree-Fock procedure is used for the 
mean- field analysis. In the ansatz-wavefunction, the 
Neel-vector is fixed in the z- and the total spin in the 

x-direction ((5) ■ (5) = 0). The pseudo-spin degrees of 
freedom of the charge/phase sector are described by an 
5 = 5 spin coherent state 
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while the spin degrees of freedom of the pair are con- 
tained in I (f> y , x) 



\4> y ,X)=e- i * y§y (cos X |A)-sin X |0». 



(11) 



\4> v ,x) i s J us t the bilayer coherent state [O where the 
global orientation of the two-spin system has been fixed. 

We list the expectation- value of a number of quantities 
with respect to the variational state 



n = 1 — (ny) = cos 2 9, 
(S x ) = nsin2xsin^ ; (S v ) = (S z ) = 0, 
(5 Z ) =nsin2xcos^; (S x ) = (S y ) = 0, 
(tia) — n cos 2 x cos 2 4> v ■ 



(G aV ) = ^/n(l-n)e-^(4>y, X \a) 



(12) 



where a — 1,0,-1,^4 labels the four spin-states. The 
role which the various parameters play can be determined 
from this list : fixes the pair-density; <p v determines the 

relative magnitude of (5) and (5), while their total mag- 
nitude is fixed by x\ ip represents the phase which orders 
in the superconducting state. 

The variational energy is given by 



E vatr .({d,t,4> v ,xh,s) = ({i,S}\H\{i,6}) , (13) 



where 



\{i,6})=U\e^;r,x)i,5. (14) 



In the mean-field analysis, it is assumed that the stag- 
gered local magnetization and the charge density are uni- 
form. The phase ipi is allowed to have a different value 
on horizontal (ip* 1 ) and vertical bonds (V' V )- We then 
arrive at the following mean-field energy 
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(sin 2 26(t 1 + t 2 - 2i 2 cos 2 x cos 2 $") 



x cos(i/' H - ip V ) - - Jf cos 4 9 sin 2 2 X 
+ i J cos 2 6(1 - 4 cos 2 x cos 2 <fi y ) - ficos 



(15) 



where N denotes the number of long bonds. Minimizing 
Eq. ([l^) , a variety of mean-field groundstates is obtained 
as a function of the various parameters. The results are 
summarized in figures ^-|^ and in table 1. We focus here 
on the case t\ > 0, for which the superconducting state 
is typically of d-wave type (ip H — i/> v = 71"). The same 
phase-diagram results for t\ — ► —t\ and i 2 — * — 1 2 , but 
with s-wave instead of d-wave phase-order. For simplic- 
ity, J, t\, t 2 and fi are expressed in units of Jp from here 
on. 

At half-filling, the physics is determined by the com- 
petition between the antifcrromagnetic and the ferro- 
magnetic spin-spin interaction. While the first promotes 
singlet-formation along the horizontal and vertical bonds, 
the second favors large local magnetic moments. J there- 
fore tunes the singlet density in the groundstate at half- 
filling. For J <C 1, the system has full Neel order with 

( n A) = hi 1(5)1 = 1- The singlet density increases lin- 
early with J up to (ua) = 1 at J = 2, where the staggered 
magnetization vanishes in a second order transition to a 
quantum paramagnet phase (Fig. 0). 
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FIG. 3. Mean field phase diagram of J versus /x and n, for 
t\ > t* (ti = OAJf, t2 — —0.1J F ). Bold lines indicate first 
order transitions. At the dotted line, transversal quantum 
spin-fluctuations destroy the antiferromagnetic order. 

FIG. 4. Mean field phase diagram of ti versus /i and n, for 
J < 2J F and t 2 < J/8 (J = J F , t 2 = -0.1J F ). 

FIG. 5. Mean field phase diagram of t 2 versus /x and n, for 
J < 2J F and h > (4j| + J 2 )/32J F (ti = 0.4J F , J = J F ). 



For densities smaller than one, the two hopping pro- 
cesses begin to play a role. Since the case of a uniform 
charge-distribution is considered and since all electrons 
are paired in the strong-coupling limit, all variational 
states with a non-integer electron-density exhibit super- 
conductivity. The superconducting order-parameter de- 
pends on the electron-density as |A| ~ yjn{l — n), see 
Eq.@. 

The value of the hopping amplitude t\ determines the 
nature of the transition from the antiferromagnetic insu- 
lator at half-filling to the singlet superconductor at lower 
densities. For small t\, this transition is first order as 
a function of /i, giving rise to a region of antiferromag- 
net/ superconductor phase separation in the t\-n phase- 
diagram (Fig. |). At ti = t\ = | + 2*1, the first order 
line splits into two second order lines. A region opens 
up in which the system has both antiferromagnetic spin 
order and superconductivity. In this antiferromagnetic 
superconductor (AFSC) phase, the electrons which carry 
the superconducting order-parameter are at the same 
time responsible for the antiferromagnetism. This state 
is most easily visualized by thinking of a small density 
of nearest-neighbor hole-pairs being doped into a half- 
filling antiferromagnet. If these hole-pairs are most mo- 
bile along the diagonals of the square lattice, where their 
movement does not disturb the AF spin order, they can 
delocalize and in that way give rise to superconductivity 
without at the same time destroying the antiferromag- 
netic order. The condition that diagonal pair-hopping 
has to dominate to get an AFSC-phase on the square lat- 
tice is reflected by the condition t\ > t\ for the present 
model. 

There are three ways in which the spin order-parameter 
in the AFSC phase is suppressed through the doping with 
hole-pairs. The simplest one corresponds with the di- 
lution of the antiferromagnet by the removal of spins. 
More interestingly, the inter-pair spin-spin interaction Jp 
scales with the pair-density squared, while the intra-pair 
spin-spin interaction J scales linearly with n. As a result, 
the ferromagnetic interaction is suppressed by a factor n 
relative to J, pushing the ratio J/Jf closer to its crit- 
ical value and reducing the magnetic moment per pair. 
Finally, the hopping process t^ frustrates the Neel order, 
provided that sign(t,2) = —sign(ti) (the other case is dis- 
cussed below). The increase of the singlet density per 
pair due to the last two processes results in a transition 
to the singlet superconductor at n = n c = gEff; ■ 

Since the *2 process amounts to a ii-type hop with 
an additional interchange of the two electrons forming 
the pair, it picks up a minus-sign when acting on a pair 
in the antisymmetric spin-singlet state. Suppose that ti 
and ti have the same sign. A singlet-pair, through the 
iVprocess, then frustrates the phase-ordering as favored 
by the ii-hop. To reduce this frustration, the singlet 
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content of the pairs is suppressed as t2 is increased, en- 
hancing the spin-ordering in the AFSC phase. Eventu- 
ally, a first order transition occurs to a ferromagnetically 
ordered triplet superconductor phase, where the singlet- 
density is reduced to zero (Fig.||). 

If t\ and ti have opposite sign, the triplet component 
is suppressed through the same process and the i2-hop 
reduces the spin-order in the AFSC phase. Note that t% 
causes a positive shift of the critical value t\ regardless 
of its sign. This reflects the fact that on-sublattice hop- 
ping must dominate in order for an AFSC phase to occur 
(ti > t\ implies t\ > ^l, where the equal sign occurs for 

\h\ = \). ' 

It should be verified that the saddle-point solution be- 
comes exact in the limit d — > oo. To reach this limit, 
the model has to be formulated in arbitrary dimension. 
The d-wave phase order then posses a problem, since it 
cannot be generalized to dimensions higher than 2. How- 
ever, the Hamiltonian of the 2d system is invariant under 
a simultaneous sign-change of t\, £2 and G^y 3 " , which 
implies that d- and s-wave order are equivalent for the 
2d model 1 . We therefore flip the sign of t\ and t^ and 
study the d — > 00 limit for the s-wave ordered state. In 
order to keep the energy finite, Jf, t\ and i 2 are scaled 
with i while taking the limit. The variation of the energy 
is found to be of order ^ for the saddle-point solution. 
Since it vanishes at large d, the mean-field groundstate 
indeed becomes an eigenstate of the system in this limit. 
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4-32(ii+t 2 ) 
1 


2n 
1 



J 

2 



TABLE 1: Mean- field results for the various phases. 



IV. TRANSVERSAL SPIN FLUCTUATIONS 

In the AFSC phase, both the U(l) phase- and the 
SU(2) spin-symmetry are spontaneously broken. As a re- 
sult, the system has two spin-wave modes and one phase 
Goldstone mode. These gapless modes dominate its long- 
wavelength physics. Since they decouple, the phase and 
spin degrees of freedom of the system may be treated 
separately at sufficiently large lengthscales. 



J This no longer holds if longer-range hopping processes are 
included. 



The physics of the phase sector is equivalent to that of 
an XY-spin model in an external magnetic field, which 
has a dynamical critical exponent z = 2. The T = 
system is therefore effectively at its upper critical dimen- 
sion d = 2 + z = 4. Because of the high effective dimen- 
sionality, phase-fluctuations only give small correction to 
the zero-temperature mean-field results for the insulator- 
superconductor transition [2^ ]. 

The long- wavelength behavior of the spin-sector in the 
AFSC phase is characterized by a critical exponent z = 1. 
Hence, at T = 0, the spin-sector lives effectively in 
three dimensions and fluctuation-effects can be signifi- 
cant. The long- wavelength spin physics is described by 
an effective non- linear sigma model p3| , p4f . This model 
contains one coupling constant, go, which is a measure 
of the quantum fluctuations in the system. At a critical 
value of go, the spin-system undergoes a quantum phase 
transition from a Neel ordered state to a quantum param- 
agnet. For the present model, go is expected to diverge at 
the mean-field transitions to the singlet superconductor 
and the paramagnetic insulator pi) . Transversal fluctu- 
ations are therefore expected to significantly reduce the 
region in the phase-diagram where AF order is stable. 

The coupling constant of the effective non-linear sigma 
model depends on the bare values of the spin stiffness and 
the perpendicular susceptibility, which are properties of 
the microscopic model. The mean-field expressions for 
these quantities can serve as an estimate for their bare 
value pi] . These expressions are derived below for the 
present model. 




FIG. 6. A spin-configuration with a twist along the 
x + y-direction. 
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We define the perpendicular susceptibility as the in- 
duced magnetization per square plaquette (containing 
four spins) by a vanishing magnetic field applied per- 
pendicular to the direction of antiferromagnetic order- 
ing. It is calculated by adding a magnetic field term to 
the mean-field energy Eq. (|l5|), 



H ( S ls) = NH cos 2 6> sin 2x sin < 



(16) 



i,S—S x ,5 y 



and subsequently minimizing the energy. This yields 



X± 



lim 

H^O H 

,AF 



2(S X ) _ 2(l-cos2x)n 



J-8t 2 (l 



X± = 



AFSC 



2-J 
J 

n(2-8t 2 )-J+8t 2 
J-8t 2 (l-«) 



(17) 



The susceptibility vanishes at the transitions to the quan- 
tum paramagnet and the singlet superconductor phase. 
It has a divergence at n = 1 — g^, which is interrupted 
by the first order transition to the triplet superconduc- 
tor phase (Fig. ||) (the line where x± diverges and the 
first-order line approach each other for small n). 

For n — 1 and Jf J ', the four spins around each 
square plaquette lock into a symmetric state. The spin- 
operators can then be replaced by 4 times the total 
spin-operator on the plaquette. The resulting Hamilto- 
nian describes an S — 2 antiferromagnet on a square 
lattice, with a spin-spin coupling J c s — A J, Such 
a system has a mean-field perpendicular susceptibility 
Xx = 1/8 J c ff = 2/J |H (where the lattice-spacing, in 
this case the distance between neighboring square pla- 
quettes, is set to one). The above result for the suscep- 
tibility indeed reduces to this expression for n = 1 and 
J< 1. 

To determine the spin stiffness, the configuration 
shown in Fig. || is considered. It has a slow twist in the 
spin order parameter along the x + ^-direction. The stiff- 
ness gives the lowest-order correction to the groundstate 
energy due to this twist p5[ . 

At each antiferromagnetic bond along the direction of 
the twist, the spins have been rotated over an angle a 8$ 
in the XZ-plane, at each ferromagnetic bond over an an- 
gle (1 — a) 8(f). This configuration is described by the 
variational state Eq. ([To|), where the spin-part is given 
by 



ju<t>s v 



1 



, X = XMF 



(18) 



with the index I labelling the bonds along the twist. 

The antiferromagnetic interaction energy of this state 
is given by (compare with last line in Eq. (|TF 



-JN cos 2 
4 



1 — 4 cos 2 x cos 2 (-a 



E AF (8(j> = 0) + -J No 2 



b 2 cos 2 9 cos 2 x : 



(19) 



while the ferromagnetic energy is simply reduced by a 
factor cos(l — a) 8<f> per twisted ferromagnetic bond. The 
phase-ordering energy also contributes to the spin stiff- 
ness. Along the twist, we have (d-wave order) 



a=l,0,-l 

- E \sm 2 20(n\e~ lLS ' i ' sv \a)(a\e^ l+1 '>^ sv \n) 

a=l,0,-l 

-isin 2 20<O|(1- \A)(A\)e i5 ' t ' sv \Q.) ~ 

\ 

(G l AV )(G l +X) - 7 sin 2 20 I 1 - -Scj? (f1|^ 2 |fi) 



/ 



(G 1 av){Gva) — 

- - sin 2 29 cos 2 v cos 2 ( -a 2 86 2 
4 A \2 

-isin 2 26»cos 2 x U - ^a 2 ^ 2 



(20) 



Taking these contributions together, the energy-increase 
due to the twist in the spin order-parameter is found to 
be 

AE = -N 84? (2 sin 2 29 [{t x + t 2 ) sin 2 X ~ ha 2 cos 2 x] 
8 

+ sin 2 2 X cos 4 6>(1 - a) 2 + 2 J a 2 cos 2 9 cos 2 x) (21) 

The distribution of the total twist over the two types of 
bonds is obtained by minimizing this energy with respect 
to a, which yields 



a 



n(l — cos 2%) 



i(l - cos2x) + J - 4t 2 (l - n) 



(22) 



For Jf much larger than J and ti (or just J for n = 1), 
the twist is entirely localized on the antiferromagnetic 
bonds, as expected. At the transition to the spin disor- 
dered phases, it is localized on the ferromagnetic bonds. 
The stiffness now follows from 



AE{a ) 



N 



(23) 



For the half- filling antiferromagnet, we obtain 
1 



„AF 



PS" = g J(2 - J) , (24) 
while for the AFSC phase the stiffness is given by 



AFSC _ n(2 - 8t 2 ) -J + 8t 2 

Ps — ~ ~ x 



8J + 16n 

x [2Jn + J 2 + 4(1 - n) (j(t x - 2t 2 ) + 2t x n + 8*§(1 - n))] 



(25) 
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Like the susceptibility, the stiffness vanishes at the tran- 
sition to a spin-disordered phase. It reduces to the 5=2, 
J off = -L j_ form for J » J F at half- filling. 

The bare coupling constant of the non-linear sigma 
model is given by go — (/OsXJ-) ^- As expected, it di- 
verges at the transitions to the singlet superconductor 
and quantum paramagnct, since both the susceptibility 
and the stiffness vanish in these phases. In order to ob- 
tain a more precise estimate of go, its value is shifted by 
a constant factor such that it agrees with the result for 
the 5=2 antiferromagnet at n = 1, Jp ^> J. The bare 
coupling for the square lattice 5=2 antiferromagnet can 
be determined from spin-wave results for the renormal- 
ized spin-wave velocity and perpendicular susceptibility, 
using the one-loop expression E4] 



.90 



1 



47r 1 + 4irx±c/%A ' 



(26) 



where A = 2^/Tr/a, with a the lattice spacing. Using the 
spin- wave results of Igarashi E3] , we obtain 



,S=2 



3.85. 



(27) 



For the half-filling antiferromagnet, the bare coupling 
constant is given by 



.90 - So " 2 



2 — J 

while we find for the AFSC phase, 



(28) 



(29) 



s=2 2y/(2n + J) (J - 8t 2 (1 ~ n)) 
g °- g ° n(2-8i 2 )- J + 8t 2 X 

1 

x — 

y/2nJ + J 2 + 4(1 - n)[J(fi - 2t 2 ) + 2t x n + 8tf(l - n)} 

The order-disorder transition at go = g c = 47r is in- 
dicated by a dotted line in the mean-field phase dia- 
grams, Fig.s [|||. It is found that transversal spin fluctua- 
tions significantly reduce the parameter-range over which 
the Neel-ordered phases are stable, without changing the 
topology of the zero-temperature phase diagram. 

At non-zero but low temperatures, the quantum non- 
linear sigma model predicts z = 1 quantum critical be- 
havior in a parameter region around the AFSC to SC 
transition line E4|. The width of this region grows as 
| <7o — 47r| ~ T~ v , with v = 0.7 the correlation-length crit- 
ical exponent of the 3d Heisenberg model. This type of 
finite temperature behavior, where temperature becomes 
the only energy-scale in the system, has been reported for 
the underdoped cuprates by a number of authors p6|-p8[ . 

Finally, we note that in the present model the super- 
conductivity onset-temperature is completely determined 
by phase-fluctuations. This is a trivial consequence of the 
fact that we focussed on the strong-pairing limit. Nev- 
ertheless, it is consistent with recent analyses of the de- 
pendence of T c on the zero-temperature phase-stiffness 



and on the number of closely-spaced layers in the super- 
conductor material, which point to a dominant role of 
finite-temperature phase-fluctuations in determining T c 

M. 



V. SO(5) SYMMETRIC POINT 

The AFSC phase has an interesting property. Let us 
consider the 50(5) superspin-vector [jllj for this model. 

Np = (^{Gav + G VA ), \s, ^(G av - G VA ^j . (30) 

The label 'P' indicates that Np is defined in the projected 
Hilbertspace, where double site-occupancy is forbidden. 
The mean-field expectation value of Np satisfies 



d(Np)* 



di 



(31) 



t2 = -- 



Hence, at the mean-field level and for this particular 
choice of t 2 , the AFSC phase can be characterized by 
an SO(5) order-parameter which has components both 
in the superconducting and in the antiferromagnetic sub- 
space, and which is rotated from the AF to the SC direc- 
tion as the hole-density is increased. As one approaches 
the tricritical point, the AFSC states with different n 
become degenerate (Fig.||) and the mean-field state be- 
comes invariant under rotations of A^p. For t 2 = — \ 
the tricritical-critical point is located at ti = t* = i, 

fJ> = V =%■ 

It should perhaps come as no surprise that we find a 
'mean-field SO(5) symmetry' for this model. The spe- 
cial lattice used here has two orbitals per unit cell, which 
seems to be one of the requirements for constructing 
an 50(5) symmetric model with short-range interactions 
ffflf . This can be understood from the fact that the min- 
imum number of sites required for the electron Hilbert- 
space in which an SO(5) representation can be con- 
structed is two (since the 7r-operators are spin-1, charge 
2 objects). Two- leg spin ladders have a natural two- 
site unit, the rung, on which the SO(5) order parameter 
can be defined. The lattice used here also has such a 
unit: the long bond. To formulate a short-range 5(9(5) 
model on the square lattice, one either has to break 
the lattice-symmetry, or to involve a certain amount of 
coarse-graining, which means that the resulting 50(5)- 
description is effective rather than microscopic. 

In the following, an exact 50(5) symmetric point is 
derived for the present model. The procedure used is 
similar to that for the 50(5) symmetric ladder (3^]. At 
the mean-field 50(5)-point, the Hamiltonian is given by 



Tt = Tio + T~ti, 



(32) 
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where 



Wo = -£ E N^-N^-JJ2 E «A> (33) 



<5 2 = 



■i S— 8 x ,5v 



^ = E E z 



5i,5i ' 5j j( 5 2 



i 6, = ±6^ 



d 2 = ±Oy 

li \ Si l S 1 ■ Si,8 2 + §i l S 1 ■ 5i.5 5 



(34) 



absorbing a d-wave staggering into the |Vi,i)-state. 

The second term, Tii, is manifestly not invariant under 
rotations of Np. After decoupling the operators on differ- 
ent bonds with respect to the order-parameters for super- 
conductivity and antiferromagnetism, this term vanishes 
and therefore the symmetry breaking does not show up 
at the mean-field level. As Eder et al. pointed out Q, 
the first and the fourth component of Np are rotated into 
each other by transforming the zero-magnetization triplet 
state |0) into the hole-pair state \V). This transforma- 
tion leaves the singlet density tia invariant. Since one 
may assume in mean-field that all components of (Np) 
vanish except the first and the fourth (spontaneous sym- 
metry breaking selects a preferred direction in the spin- 
and phase sector) the decoupled mean-field Hamiltonian 
is invariant under this transformation. This implies that 
the d — > oo 50 (5)-symmetry is not only present in the 
zero-temperature groundstate, but also at finite tempera- 
tures, where higher energy-levels are thermally occupied. 

As a first step towards an 50 (5)-symmetric Hamilto- 
nian, TL\ is subtracted from TL. This introduces second- 
and third-neighbor spin-spin interactions into the model. 

The second term in Tlo is an 50(5)-invariant (this is 
discussed below). The first term is invariant under ro- 
tations of Np , but this does not imply that it is 5*0(5) 
symmetric. There is no representation of the 50 (5) alge- 
bra on the projected Hilbertspace under which Np trans- 
forms as a vector. The rotation-symmetry is therefore 
broken at the quantum level. In a recent article p5[ , 
Zhang et al. show that mean-field 50(5) symmetry al- 
ways remains when a projection to the lower Hubbard 
band is applied to a system with full 5*0(5) symmetry. 
Here we work backwards: mean-field 50(5)-symmetry 
being established, we deduce a model with full SO (5) 
symmetry by lifting the constraint of no double site- 
occupancy. 

The basis of the single-bond Hilbert-space is extended 
to include the doubly occupied state \D). It now consists 
of one 50(5) singlet (A) and one 50(5) quintet (spin- 
triplet, D and V). The details of this representation of 
the 50(5) algebra are briefly discussed in appendix |A|. 
We introduce an on-site repulsion U J^is n DiS- The gen- 
eral Hamiltonian on the unprojected Hilbert space is 
given by 



* = E E 



4(* 1 +i 2 )E^ 1 



+(ti - t 2 )A i;Sl ■ A ii4a 



(35) 



E 

i,S=5^,S y 



Jn A S + (T - U + V) n D ~ ~ T) 



i.S 



i,6 



where A = (ReA,ImA), 7?" = (Re7r Q , Im7r Q ) (see ap- 
pendix |A|) . The value of U has to be fine-tuned in order 
to obtain 50(5) symmetry on a single bond. The re- 
sulting constraint is fj, = i as before, but in addition 
f7 = ~. Note that this is more restrictive than the local 
constraint for the ladder model |nj , which leaves two free 
parameters. Since we only consider states of paired elec- 
trons, this model has fewer local 50(5) invariants than 
the ladder. 

To establish 50(5) symmetry of the inter-pair interac- 
tions, one now has to take t\ = — t 2 = \- After subtrac- 
tion of Tii, this yields the Hamiltonian 



H 



SO(5] 



= -7 E NrN m -jJ2' 



■'A ■ 



(36) 



<l,m> 



where / and m run over the square lattice spanned by 
the long bonds (dotted lines in Fig.|lj) . The unprojected 
50(5) superspin N is given by Eq. (|Al|). The local 50(5) 
invariant tia is related to the length of the superspin 
through iV 2 = 1 + An a- 

Note that Eq.(|36]) is not the most general 50(5)- 
symmetric Hamiltonian which could be formulated. In 
principle, there can be an additional term of the form 



a E E^s = a E 



</,m> a<b 



<l,m> 



+5i i 5 1 ■ Si, 



Qi,S 2 



(37) 



which is also an 50(5) invariant. The charge-charge in- 
teraction Qi,8 1 Qi,s 2 was omitted from the present analy- 
sis and a term of this form therefore does not appear at 
the symmetric point. 

The projected 50(5) (p50(5)) symmetry at U — > oo 
evolves from the true 50(5) symmetric point at fine- 
tuned U in the following way. Let us assume we have 
ti = -t 2 = t and U = i + U. 



H 



E ( 2< A; • A T 

<l,m> 



1 



-Si ■ S r7 



+ E ( Un D ~ Jn 



(38) 



The superspin has no preferred global direction for U — 
0, t = |. Since the AF groundstate does not have a \D) 
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component, while the SC does, a small positive U will flop 
the superspin to the AF direction. The energy-difference 
between the AF and the SC state can be compensated 
by an increase in t. For U — > oo, this procedure shifts 
the superspin- flop point from t = g to t = -j, with the 
SC groundstate now having {no) = 0. The shift in t is 
accounted for by the different relative normalization of 

A and 5 in the definitions of Np and N. 



VI. COLLECTIVE MODES 

In the above, it was shown that the inter-sublattice 
hopping ti , which couples the spin- and charge-dynamics 
in our model, plays a crucial role in establishing the 
50(5)-symmctry. This symmetry only emerges at the 
mean- field level for a specific value of t 2 - To further 
investigate the role of this hopping process, the collec- 
tive modes in the antifcrromagnetic and spin-disordered 
phases are analyzed. For t 2 = 0, it is found that a decou- 
pled spin/charge perspective suffices to understand these 
modes, as one would expect. In this case, the supercon- 
ductivity does not affect the collective spin-modes of the 
system. 

This changes for non-zero t 2 . Although the dispersion- 
relations do not change qualitatively, the interpretation 
of the modes does. Most strikingly, the gapped spin- 
magnon mode of the Singlet SC phase acquires a 7r-mode 
component. This mode softens at the transition to the 
AFSC phase and becomes a pure, acoustic 7r-mode as the 
system is tuned towards p50(5)-symmetry. 

The mode-spectrum for systems with p50(5)- 
symmetry was analyzed by Zhang et al. fl5|| . We repro- 
duce their results for the present model and investigate 
the influence of further 5*0(5) symmetry-breaking terms. 



A. Random phase approximation 

The collective modes of the system are studied in the 
random phase approximation (RPA) |Q. They are ob- 
tained from the equations of motion of the operators G a p, 
which are given by 



(39) 



The commutator in this expression contains products of 
operators on different bonds (i,S). These products are 
decoupled in a mean- field fashion, yielding a set of cou- 
pled linear differential equations. After a transformation 
to frequency- and momentum-space, it takes the form 



uiG a(} {k,ui) 



j2M:'f(k)G a , 0l (k, t 

a'B' 



(40) 



while its eigenvectors give the operators which generate 
these modes. 

There is a problem with the above decoupling in the 
spin-ordered phases. As was discussed in section [V, the 
low-energy fluctuations of the spin-system behave differ- 
ently at large Jf and near the spin-disordering transition. 
For Jp <C 1, the spins around each square plaquette are 
locked into a symmetric state, forming one spin-2 ob- 
ject, and the low-energy deformations of the spin-state 
are localized on the antiferromagnetic bonds. The above 
decoupling, which cuts across the ferromagnetic bonds, 
then becomes very poor. Near the spin-disordering tran- 
sition, the spins are rigidly coupled along the antiferro- 
magnetic bonds and the low-energy transversal fluctua- 
tions are localized on the ferromagnetic bonds. In this 
case, the decoupling works well. 

The cross-over to spin-2 behavior at large Jp is driven 
by the 7i\ spin-spin interaction term, Eq. (|34|) . To avoid 
it, we calculate the mode-spectrum of the Hamiltonian 
from which this term has been subtracted. The result- 
ing dynamical matrices are listed in appendix Sub- 
tracting Tii makes no difference for the spin-disordered 
phases, but does change the results in the AF and AFSC 
phase (we discuss in what way). The TL\ term breaks 
50(5) symmetry, though retaining it at the mean-field 
level. As a result, the model which we study in RPA has 
a projected 50(5) symmetry at the tricritical point with 
fine-tuned t 2 . 




The dispersion-relations of the collective modes are ob- 
tained from the eigenvalues of the dynamical matrix M, 



FIG. 7. Sketch of the mode-spectrum in the antiferromag- 
netic and spin-disordered phases. The dashed lines are pairing 
modes when gapped (insulating phases) and phase Goldstone 
modes when acoustic (superconducting phases). 



B. mode-spectrum 

We briefly discuss the mode-spectrum of the antiferro- 
magnetic and spin-disordered phases. These results are 
summarized in Fig. [?] . 

The quantum paramagnet, for which both spin- and 
gauge-symmetry are unbroken, has no Goldstone modes. 
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Its spectrum consists of a gapped three-fold degenerate 
spin-1 magnon mode and a pairing mode which is also 
gapped. The pairing gap closes at the transition to the 
singlet superconductor. Precisely at the transition, this 
mode has a quadratic dispersion. As the hole-pair den- 
sity is increased from zero, it acquires a finite velocity, 
becoming the phase Goldstone mode of the supercon- 
ducting state. This is the behavior expected at a dilute 
boson transition p2[ , of which this is a particular ex- 
ample (where the hole-pairs are the dilute bosons). The 
three-fold degenerate mode of the quantum paramagnet 
remains gapped through the transition to the supercon- 
ductor. 

The insulating antiferromagnet has a two-fold degen- 
erate spin-wave mode. In addition, it has a gapped mode 
related to spin-amplitude modulations and a gapped 
pairing mode. The spin-amplitude mode becomes de- 
generate with the acoustic spin-wave modes at the tran- 
sition to the quantum paramagnet, where they turn into 
the spin-1 magnon triplet. This transition is in the same 
class as the spin-disordering transition of the Heisenberg 
bilayer model (3^]. The pairing gap closes at the tran- 
sition to the antiferromagnetic superconductor, where it 
becomes the acoustic phase Goldstone mode at finite dop- 
ings. The spin-amplitude mode remains gapped through 
the insulator to superconductor transition, but becomes 
degenerate with the acoustic spin-wave modes at the sub- 
sequent transition to the singlet superconductor. At this 
transition line, the system therefore has four acoustic 
modes, of which three are degenerate. 



C. 7r- modes 

The mode-spectrum as outlined above can be under- 
stood entirely from a decoupled perspective of the spin- 
and the phase-sector. At the transition line from insu- 
lator to superconductor, the pairing mode softens and 
continuously acquires a finite velocity. At the transition 
from a phase with AF order to a spin-disordered phase, 
the spin-amplitude mode becomes degenerate with the 
acoustic spin-wave modes. These two effects combined 
yield the described behavior, and in particular the oc- 
currence of four acoustic modes at the AFSC to Singlet 
SC transition. This decoupled perspective is correct for 
*2 = 0. In this case, the gapped modes in the Singlet SC 
state are indeed spin-1 magnons, as they are in the quan- 
tum paramagnet, while the gapped mode in the AFSC is 
indeed a pure spin-amplitude mode, as it is for the AF 
insulator. 

The t 2 -process, however, provides a coupling between 
the movement of the hole pairs and the dynamics of 
the spin-system. This coupling changes the nature of 
the gapped modes in the superconducting phases with 
respect to those in their insulating parent-phase. In 
the AFSC phase, the spin-amplitude mode is mixed 



with fluctuations between the hole-pair and the zero- 
magnetization electron-pair states. In the Singlet SC 
phase, 7r-modes (hole-pair to triplet fluctuations) are 
mixed into the spin-1 magnons. 

Close to the transition from the Singlet to the anti- 
ferromagnetic SC, the gapped mode becomes degenerate 
with the acoustic spin-wave modes and the ^-process be- 
gins to affect the low-energy physics of the system. From 
the analogy with the spin-disordering transition in the 
insulating phase, one would expect to find a three-fold 
degenerate acoustic magnon-mode at this transition (in 
addition to the phase mode). Instead, the RPA-analysis 
yields an eigenvector 



4t 2 VJ^2{G v-Gvo) 

+(J - 8* 2 )(1 - 2ti + 2* 2 ) (Gqa - Gao) , (41) 

which has both a magnon and a 7r-mode component. 
Note that this result does not change if the more natural 
spin-spin interactions, with Tii, are used, since H.i does 
not affect the collective modes in the singlet SC phase. 

For J — > 2, the spin-disordering line in the super- 
conducting phase approaches half-filling. In this case, 
Eq.([4l|) becomes a pure magnon-mode, which is the re- 
sult expected for the insulating phase. The same eigen- 
vector is found for * 2 = 0, which implies that the tran- 
sition in the superconducting phase is, for that case, in- 
deed of the same type as for the insulators. The 7r-modcs 
are mixed in for finite *2. We find a pure 7r-mode for 
* 2 = J/8 and *i — * 2 = 1/2. The first condition is satis- 
fied at the point where the AFSC, triplet SC and singlet 
SC meet, at n — > (see FigJ|). This is related to the 
fact that the n = 0-state and the triplet SC, which are 
related by a 7r-rotation, become at that point degener- 
ate in energy. The second condition is of more interest: 
it is fulfilled if the singlet-pair hopping-process and the 
Neel-moment interaction enter the Hamiltonian in the 
projected 50(5)-symmetric form Np' Sl ■ Np S2 . This is of 
course the case at the pS'0(5)-point, which implies that 
the Singlet SC phase at this point has the mode-content 
expected from S'0(5)-theory: 4 acoustic modes, of which 
one is a phase-mode and three are 7r-modes. It is shown 
in Ref. [3 that this is generally the case for systems 
with a projected 5*0(5) symmetry in the SC phase. The 
symmetry breaking due to the projection onto the lower 
Hubbard band shows up in the RPA mode-spectrum by 
a different velocity for the phase- and the 7r-modes: 



v 



SC,pSO(5) 
phase 

SC,pSO(5) 



J 2 



4V2 



(42) 



The two modes become degenerate at J=6/5 (n=4/5), 
but this point does not seem to have any special signifi- 
cance, 
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The condition t\ — t% = 1/2 can also be satisfied at the 
AFSC to Singlet SC transition away from the point with 
mean-field 5*0(5) symmetry. In this case, there are addi- 
tional terms which break the mean-field SO (5) symme- 
try, since they tune the system away from the tricritical 
point, but which do not affect the RPA-modes. 



D. Projected SO(5) symmetry 

In Ref. |HJ, the mode-spectrum at the pS'0(5)-point 
was studied for a general direction of the superspin. Two 
striking results were obtained. In the first place, the 
system has a two-fold degenerate acoustic mode, whose 
velocity is independent of doping (i.e., independent of the 
direction of the 5'0(5)-order parameter). Secondly, the 
phase Goldstone mode is not acoustic, but gapless with 
a quadratic dispersion. It was argued that this last effect 
is caused by the infinite compressibility of the system at 
the pSO(5) point, where d(n)/d/j, diverges. 

Both results are reproduced in our model. Since the 
second result is related to the infinite compressibility 
rather than to S f O(5)-symmetry, it always occurs at the 
tricritical point, also if we tune away from t\ = — t% = j 
while keeping t\ = t*,fi = fi c . By the same argument, 
this effect will not disappear if the Hi -term is added to 
the Hamiltonian, since this does not affect the mean-field 
phase-diagram. 

The first result is very sensitive to perturbations. As 
soon as the tricritical point is tuned away from ti — 
—ti = \, the velocity of the acoustic modes becomes n- 
dependent. Also the addition of 7i\ to the Hamiltonian 
destroys this effect. This can be seen by calculating the 
spin-wave velocity from c = yj p s /\±, using the results 
obtained in section [V, and evaluating it at the mean- 



field SO(5)-point. This yields 



„mfso( 5 )_ / (2 + J)(J + 2-2n))(J+l-n) 



8(J+2n) 



(43) 



which has an n-dependence. For the model without Hi, 
the stiffness is given by (see Eq.(E3p3)) 



P. 



H-Hi 



N6(f> 2 
n(l — n)(ti 



AE{a = 0) 



(44) 



t 2 )(l-cos2 X ) + T (l-cos 2 2 X )], 



where n and \ have the mean-field values listed in ta- 
ble 1. The susceptibility is obtained by multiplying the 
J^-term in the mean-field energy Eq. (|l5|) with a factor 
cos 2 <t> v , adding the magnetic field term Eq.([l6|), and min- 
imizing with respect to </> y . We obtain 



2rc(l - cos2x) 



At the pSO(5) point, this stiffness and susceptibility re- 
produce the RPA result v s = (2 + J) /(Ay/2), which is 
independent of doping. Both the phase-ordering and the 
spin-ordering energy contribute to the spin-wave veloc- 
ity in the AFSC phase. As one approaches the pSO(5) 
point, the doping-dependence of the contribution of the 
spin-ordering energy is precisely compensated by the op- 
posite doping-dependence of the phase-ordering contri- 
bution. Note that the model with Tti yields the same v s 
at the transition to the Singlet SC, where n = (2 + J)/4, 
see Eq. (fi"3|). This demonstrates the insensitivity of the 
RPA mode-spectrum in the Singlet SC phase to the spe- 
cific form of the spin-spin interactions. 

At the pSO(5) point of our model, the acoustic modes 
are no longer pure spin- wave, but a combination of spin- 
wave and 7r-mode. Their eigenvector is given by 



J - 8i 2 (l - n) + n(l - cos 2 X ) 



(45) 



VT^(Giv + G v _i) + y/h~^h~ c (G w + G _i) , (46) 

which starts out as a spin-wave at half-filling, but crosses 
over to a 7r-mode at n = n c . This agrees with Zhang et 
aVs interpretation of the doping-indcpcndcnt velocity in 
terms of the projected ,90(5) symmetry [jisj] . 



VII. SUMMARY 

We have introduced a strong coupling model for spin- 
ordering and superconductivity. The microscopic build- 
ing blocks of this model are nearest-neighbor electron 
pairs. The spatial structure of these pairs gives rise to 
d-wave superconductivity. At the same time, it allows 
the pairs to have a non-zero uniform or staggered mag- 
netic moment. In order to avoid problems related to 
dimer-type spatial correlations between the pairs, the 
model is formulated on a 1/5-depleted lattice. A rich 
mean-field phase-diagram is obtained, exhibiting in par- 
ticular a phase which is at the same time an antiferro- 
magnet and a superconductor. The second order lines 
separating this phase from the antiferromagnetic insu- 
lator and the spin-disordered superconductor end at a 
tricritical point, where the antifcrromagnet to supercon- 
ductor phase-transition becomes first order. By map- 
ping the spin-sector in the antiferromagnetic phases onto 
a non-linear sigma model, the main corrections to the 
mean-field phase-diagram have been obtained. 

For a specific value of one of the model parameters, a 
mean-field ^0(5) symmetry between the antiferromag- 
netic and superconducting order-parameter appears to 
be realized at the tricritical point. It turns out that 
the model still contains spatial gradient terms which 
break SO(5) symmetry. These can be removed by mod- 
ifying the spin-spin interactions. The remaining 5*0(5) 
symmetry-breaking is then a pure quantum-effect, being 
related to the operator-algebra rather than the Hamilto- 
nian. It is shown that true £0(5) symmetry can be real- 
ized for this model by allowing double site-occupancy and 
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fine-tuning the Hubbard U. The approximate symmetry 
at large U is therefore a projected SO (5) symmetry of 
the kind discussed in Ref. 

We investigated the mode-spectrum using the random 
phase approximation. It is found that the inter-sublattice 
hopping process gives rise to the appearance of a n- 
component in the gapped modes of the Singlet SC phase. 
Approaching the point with projected SO(5) symmetry 
from the Singlet SC phase, a three-fold degenerate acous- 
tic 7r-mode is found as well as an acoustic phase-mode. 
The RPA mode-spectrum then has the properties ex- 
pected for an 50(5)-symmetric system in the pure super- 
conducting phase, apart from the fact that the 7r-modes 
and the phase-mode have different velocities. 

As reported in Ref. fl5|] , the system at the pSO(5) 
point has a gapless phase-mode with a quadratic dis- 
persion, as well as a two-fold degenerate acoustic mode 
whose velocity is independent of doping. This acoustic 
mode crosses over from a pure spin-wave at half-filling 
to a 7r-mode at the transition to the Singlet SC phase. 
Wc find that the quadratic phase-mode is a property of 
the tricritical point rather than of the projected SO (5) 
symmetry. The doping-independent velocity, however, is 
a strong signature of projected SO (5) symmetry, which 
can be destroyed even by additional symmetry-breaking 
terms that leave the mean-field 5*0(5) symmetry intact. 
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(A4) 



where the 7r-operators read iv a 



c\a a a v c\, with a the 



Pauli matrices jjl). Projecting onto the paired-electron 
states, we obtain 



A - 



— (Gdi — Gd -i + G\v — G-xv) 



■ = — (Gdi + Gd-i — Gw — G-iv) , 



7r| = -j= (Gdo + Gov) ■ 

The charge-operator is given by 
Q = tid — ny . 



(A5) 



(A6) 



It can be checked that N indeed transforms as a vector 
under this S f O(5)-algebra: 



[L ab , N c ] = i (S ac N b - SbcNa) 
and furthermore that 

[N a ,N b ] =iL ab . 



(A7) 



(A8) 



APPENDIX A: THE SO (5) ALGEBRA. 



APPENDIX B: DYNAMICAL MATRICES 



A short overview is given of the representation of the 
50(5)-algebra for this model. 

In the unprojected Hilbcrtspacc, a representation of 
the 5*0(5) algebra can be defined which transforms the 
superspin N as a vector. The superspin is given by 



where 



N = (ReA, §, ImA) , 



At =V2(G DA ~G AV ) 



(Al) 



(A2) 



and ReA = \ (A^ + A), ImA = i (At - A) . The gen- 
erators of the 50(5)-algebra satisfy the commutation re- 
lation 

[L a b, L c d] — 1 (SacLbd + SbdLac ~ 5adL bc — 5 0c L a d) , 

(A3) 

where the indices take the values 1 through 5. The L ab 
are anti-symmetric under an interchange of a and b. They 
are given by pit] 



A staggering factor for the antiferromagnetic spin- and 
the d-wave phase-order has been absorbed into the op- 
erators G a f3. After subtraction of Tii, the Hamiltonian 
takes the form of a model on the square lattice, where the 
operators G a p act on the states on the lattice sites. The 
Singlet dSC, AF dSC, quantum paramagnet and AF in- 
sulator mean-field states are all uniform in terms of these 
operators. It is therefore not necessary to introduce a 
multi-sublattice structure. The different modes in terms 
of the real (non-staggered) operators are simply related 
to the ones obtained here by a shift in fc-space. 

The operators G a ^ separate into three sets. Each op- 
erator couples only to operators in the same set through 
its RPA equation of motion. One set is formed by the 
raising operators {Gi V , Gia, G w , G V -i, Ga-i, G -i}, 
another by the lowering operators, which are related to 
this set by hermitian conjugation. The third set contains 
the operators which act only on the zero-magnetization 
states: {Gav, G v, G ao , Gva, Gvo, Gqa, n Al n 0l ny}- 

The dynamical matrix of the raising operators has the 
form 
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Me 



A T B T 
-B T -A T 



where A and B are the 3x3 matrices 



A = 2s 2 ec x {Y^ak - A t ) 
V 2E t s 2 0S x (l - 7 fe ) 



-2s 2e c x (A t - 2f) 



(Bl) 



£ t 2s 2 0S x 



5 



^s 2e c x 7 fc 
cec x 7fc 







(B2) 



(B3) 



2 ces 2x7fc 

where we have used the notation 

s x = sin x ; c x = cos x , 
St = t\ + t% ; Aj = t\ — t2 ■ 
1 



(B4) 



7fc = -^(cos(fc • ei) + cos(fc • e 2 )) . 



with e\ and e 2 the basis-vectors of the square lattice 
spanned by the long bonds, Fig. [j]. The angles \ and 6 
are the ones appearing in the mean field energy Eq.(15). 
In the insulating phases, 6 vanishes, while \ is equal to 
zero in the spin-disordered phases. 

The dynamical matrix of the lowering modes is the 
same as Mr, apart from a minus sign. The last set has 
a dynamical matrix 



-2Ats 2fl c x (l- 7fc ) 

E = I 2S t s 2 es x (l-7fc) 

2A t s 2e c x (l - 7 fc ) -2S t s 2e s x (l - 7fe ) 







F = 



-2A t s 2 gc x 


2Ats 2( 9C x 



c^s 2x 
2E t s 2e s x -c 2 g s 2x 



3 x 

-2£ t s 2 0S x 



(B8) 



(B9) 







For the Singlet SC phase (\ = 0), the operators 
G a fj with a and (3 referring to triplet states decou- 
ple from the equations of motion, since there is no 
longer a triplet component in the mean-field groundstate. 
This leaves the raising set {Giy, G\a, Gv-i, Ga-i} 
and the hermitian conjugate lowering set. The set of 
zero-magnetization operators splits into the transver- 
sal set {Gov, Gqa, Gvo, Gao} and the longitudinal set 
{Gav, Gva, nA — fiv}. The dynamical matrices of the 
first three sets contain the three-fold degenerate n- and 
spin-1 magnon modes. The fourth set contains the phase 
and the pairing mode. 

In the quantum paramagnet phase (x = = 0), the 
first three sets further simplify. Since the 7r-operators do 
not refer to the spin-singlet groundstate, only the spin- 
l/charge-0 operators are left in these sets. The phase- 
mode tia — ny disappears from the last set. 



M 



C D T E T 
~D T -C -E T 
F -F 



which consists of the 3x3 matrices 



D = 



C n = -4A t (c x c^ 



Cl2 
Cl 3 
G21 

C22 
C23 
C31 

C32 
C33 



4hk 



2S t c e s 2x 7fc - c s 2x 

-4E*((44 - 4Yt h - 


2E f s 2e s x - -s 2 9S x 



1 



; s 2ec x7 fc 
-CgC 2x 7fc + J 



(B5) 



2A t Cgs 2x 7fe - c 2 g S2 X 
-2s 2 gs x (A t 7fe - S t ) 



^4 



(B6) 



2 S 20S x Jk 



-2s 2e c x (A t - S t 7fc) 



2AtS 2e c x - 2 s 2ec x 7fe 



c e c 2 X 7fe 
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